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Abstract 


We  present  an  approximation  algorithm  for  the  traveling  salesman  problem 
when  the  vehicle  is  constrained  to  move  forward  along  paths  with  bounded  cur¬ 
vatures  (Dubins’  vehicle).  A  deterministic  algorithm  returns  in  time  O(n^)  a  tour 
within  ^^1  +  max  |  logn^  of  the  optimum  tour,  where  n  is  the  num¬ 

ber  of  points  to  visit,  p  is  the  minimum  turn  radius  and  Dmin  is  the  minimum 
Euclidean  distance  between  any  two  points.  A  randomized  version  returns  a  tour 
with  an  expected  approximation  ratio  of  (  ( 1  +  )  log  n  ) .  This  very  simple  al- 

y  y  ^min  J  J 

gorithm  reduces  the  Dubins  traveling  salesman  problem  to  an  asymmetric  traveling 
salesman  problem  on  a  directed  graph. 


1  Introduction 

In  the  Traveling  Salesman  Problem  (TSP),  we  are  given  n  nodes,  and  for  each  pair  (i,  j) 
of  distinct  nodes,  a  distance  dij.  We  want  to  hnd  a  closed  path  that  visits  each  node 
exactly  once  and  incurs  the  least  cost,  which  is  the  sum  of  the  distances  along  the  path. 
The  distances  need  not  be  symmetric,  i.e.  we  can  have  two  nodes  with  ^  dji, 
in  which  case  the  problem  is  called  the  asymmetric  traveling  salesman  problem  (ATSP). 
In  the  metric  TSP,  the  distances  satisfy  the  triangle  inequality.  A  subcase  of  the  metric 
TSP  is  the  (planar)  Euclidean  TSP  (ETSP),  where  the  nodes  he  in  and  the  distance 
is  the  usual  Euclidean  distance. 

Computing  good  TSP  tours  efficiently  is  of  interest  in  the  area  of  aerial  surveillance. 
As  we  are  increasingly  interested  in  developing  autonomous  vehicles,  the  question  of  how 
these  vehicles  should  behave  usually  leads  to  optimizing  a  given  objective  function,  for 
example  minimizing  the  distance  traveled  when  the  task  is  to  explore  a  set  of  locations. 
An  important  difficulty  arises,  however,  when  the  problem  involves  planes,  underwater 
vehicles,  cars  and  other  vehicles  with  significant  dynamics:  the  paths  obtained  from 
algorithms  solving  the  Euclidean  TSP  are  infeasible.  Kinodynamic  planning  refers  to 
the  path  planning  problem  when  the  kinematic  constraints  of  the  vehicle  are  taken  into 
account.  The  methods  developed  in  this  held  aim  at  Ending  a  trajectory  from  an  initial 
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position  and  configuration  to  a  final  position  and  configuration,  usually  while  avoiding 
potential  obstacles.  In  this  paper,  we  study  a  different  problem.  We  want  to  optimize 
trajectories  visiting  a  specihed  set  of  points,  but  the  conhguration  of  the  vehicle  at  these 
points  is  free  as  long  as  the  kinematic  constraints  are  satished. 

We  focus  on  algorithms  computing  tours  for  the  Dubins’  vehicle,  a  problem  which  was 
recently  considered  in  [10] .  The  Dubins’  model  [3]  is  a  simple  but  efficient  way  to  handle 
the  dynamic  characteristics  of  hxed-wing  aircrafts.  It  gives  a  complete  characterization 
of  the  optimal  paths  between  two  conhgurations  for  a  vehicle  with  limited  turning  radius 
moving  in  a  plane  at  constant  speed. 

The  complexity  issue  of  our  algorithms  is  an  important  part  of  the  analysis,  since 
the  path  planner  is  usually  only  one  component  of  a  general  scheduling  system,  and 
therefore  its  task  should  not  become  prohibitively  time-consuming.  The  TSP  is  NP-hard, 
even  in  the  Euclidean  case  [8],  and  therefore  we  are  interested  in  efficient  approximation 
algorithms.  More  precisely,  an  a- approximation  algorithm  (a  >  1)  for  a  minimization 
problem  with  optimum  OPT  is  an  algorithm  that  produces  in  polynomial  time  a  feasible 
solution  whose  value  Z  is  within  a  factor  a  of  the  optimum,  i.e.  such  that 

OPT  <Z<a  OPT. 

It  is  known  that  in  the  most  general  case,  there  can  be  no  approximation  algorithm  for 
the  TSP  unless  P=NP.  But  if  the  distances  satisfy  the  triangle  inequality,  Christohdes’ 
algorithm  [2]  gives  a  3/2-approximation  for  the  symmetric  TSP,  and  there  is  a  (logn)- 
approximation  for  the  ATSP  [4]  (where  n  is  the  number  of  points,  and  log  denotes 
the  logarithm  of  base  2.  No  constant  factor  approximation  for  the  ATSP  is  currently 
known).  For  the  ETSP,  Arora  [1]  gave  a  polynomial-time  approximation  scheme  that 
can  approximate  the  optimal  tour  within  (1  -|-  e)  for  any  e  >  0.  The  analysis  of  our 
algorithm  will  show  that  we  have  a  Oi  ^  approximation  for  the  Dubins  TSP,  where 

^min 

Drain  is  the  minimum  Euclidean  distance  between  any  two  points  in  the  set. 

This  paper  is  organized  as  follows:  in  section  2,  we  recall  some  facts  about  point-to- 
point  Dubins’  paths.  Then  in  section  3,  we  present  a  simple  approximation  algorithm  for 
the  Dubins’  traveling  salesman  problem  (DTSP)  and  compare  it  to  the  approach  taken 
in  [10].  Section  4  gives  a  randomized  version  of  the  algorithm  that  achieves  a  slightly 
better  approximation  factor.  Section  5  presents  computational  results,  and  hnally  we 
conclude  in  section  6  on  the  tightness  of  our  analysis  and  possible  improvements. 

2  Point-to-Point  Dubins’  Paths 

We  consider  a  point  vehicle  moving  at  unit  speed  (without  loss  of  generality)  in  the 
plane,  with  a  constraint  on  its  maximal  turning  rate.  More  formally,  given  a  continuously 
differentiable  path  P  :  /  — *>  parametrized  by  arc  length  s  G  /,  the  average  curvature 
of  P  in  the  interval  [si,S2]  C  /  is  dehned  by  ||P'(si)  —  P'(s2)||2/|si  —  S2I.  We  require 
that  the  average  curvature  of  the  vehicle’s  path  be  at  most  p  in  every  interval.  Denote 
the  configuration  of  the  vehicle  by  {X,6),  where  X  is  the  location  of  the  vehicle  in  the 
plane  and  6  G  (— tt,  tt]  is  its  heading,  i.e.  the  angle  that  the  velocity  vector  makes  with 
the  x-axis.  The  goal  is  to  design  an  efficient  approximation  algorithm  which,  given  a 
set  of  point  locations  in  a  bounded  square  (in  an  obstacle- free  environment),  returns  a 

^We  say  f{n)  =  0{g{n))  if  there  exists  c  >  0  such  that  f{n)  <  cg{n)  for  all  n,  and  f{n)  =  fl{g{n))  if 
there  exists  c  >  0  such  that  /(n)  >  cg{n)  for  all  n. 


permutation  of  the  points  specifying  the  order  of  the  visits,  as  well  as  headings  for  the 
vehicle  at  each  point. 

Dubins  [3]  characterized  curvature  constrained  shortest  paths  between  an  initial  and 
a  hnal  conhguration.  Let  P  be  a  feasible  path.  We  call  a  nonempty  subpath  of  P  a 
C-segment  (resp.  S'-segment)  if  it  is  a  circular  arc  of  radius  p  (resp.  a  straight  line 
segment).  We  paraphrase  the  following  result  from  Dubins: 

Theorem  1  ([3]).  An  optimal  path  between  any  two  configurations  is  of  type  CCC  or 
CSC,  or  a  subpath  of  a  path  of  either  of  these  two  types.  Moreover,  to  be  optimal,  a  CCC 
path  must  have  its  middle  arc  of  length  greater  than  np. 


In  the  following,  we  will  refer  to  these  minimal-length  paths  as  Dubins’  paths.  When  a 
subpath  is  a  C-segment,  it  can  be  a  left  or  a  right  hand  turn:  denote  these  two  types  of  C- 
segments  by  L  and  R  respectively.  Then  we  see  from  theorem  1  that  to  hnd  the  minimum 
length  path  between  an  initial  and  a  hnal  conhguration,  it  is  enough  to  hnd  the  minimum 
length  path  among  six  paths,  namely  among  {LSL,  RSR,  RSL,  LSR,  RLR,  LRL}.  Each 
of  these  paths  can  be  explicitly  computed  (see  for  instance  [11])  and  therefore  hnding  the 
optimum  path  and  its  length  between  any  two  conhgurations  can  be  done  in  constant 
time. 

Note  that,  in  our  case,  the  conhguration  of  the  vehicle  at  each  point  is  not  completely 
specihed.  Only  the  position  of  the  points  is  known,  the  headings  of  the  vehicle  at  these 
points  must  be  found.  Since  we  do  not  expect  to  obtain  the  exact  optimum  headings, 
we  will  need  a  lemma  based  on  a  result  due  to  Jacobs  and  Canny  [6]  on  the  diherence  in 
length  between  two  paths  when  only  the  initial  and  terminal  headings  are  diherent.  In 
general,  if  we  consider  a  path  of  Dubins’  length  pij  between  two  conhgurations  (Xj,  6i) 
and  {Xj,6j)  and  make  an  error  up  to  J  G  (— 7r,7r]  on  the  initial  and  hnal  headings,  we 
can  derive  a  multiplicative  bound  on  the  perturbed  path  pij  as  follows. 

Lemma  2.  Let  dij  be  the  Euclidean  distance  between  Xj  and  Xj.  We  have: 


pij  <  1  +  2pmax 
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See  [6]  for  a  proof  of  this  lemma.  Let  us  mention  that  the  two  terms  in  the  max  on  the 
right-hand  side  correspond  to  the  cases  where  the  initial  Dubins’  path  is  a  CSC  path  with 
opposite  initial  and  hnal  turning  directions  and  a  CCC  path  respectively.  Perturbations 
of  a  CSC  path  with  identical  initial  and  hnal  turning  directions  are  dominated  by  the 
hrst  term. 


3  An  Approximation  Algorithm  for  the  DTSP 

Savla  et  al.  [10]  gave  an  algorithm  for  computing  Dubins’  TSP  tours,  called  the  “alter¬ 
nating  algorithm”.  It  works  as  follows:  given  a  set  of  n  points,  the  optimal  Euclidean 
TSP  tour  is  computed,  and  the  order  of  visits  for  the  ETSP  is  used  for  the  Dubins’  tour. 
It  is  then  necessary  to  obtain  a  feasible  path  through  these  ordered  points.  Following 
the  edges  of  the  tour,  all  odd-numbered  edges  are  retained  (i.e.  the  subpath  is  a  straight 
line)  as  well  as  the  corresponding  headings,  and  the  even-numbered  edges  are  replaced 
with  Dubins’  paths  (see  Fig.  1). 

The  alternating  algorithm  is  not  an  approximation  algorithm  for  the  Dubins’  TSP  in 
the  sense  of  the  definition  given  in  the  introduction.  It  could  easily  be  made  a  polynomial 


Alternating  Algorithm 


Figure  1:  DTSP  tour  obtained  from  the  alternating  algorithm  for  a  set  of  10  points. 


time  algorithm  by  using  Arora’s  (1  +  e)-approximation  of  the  ETSP  instead  of  an  exact 
algorithm.  But  a  more  critical  issue  is  that  there  could  be  little  relationship  between 
the  Euclidean  and  Dubins’  metrics,  especially  when  the  Euclidean  distances  are  small 
with  respect  to  the  turning  radius,  which  is  the  case  we  are  mainly  interested  in  for  the 
Dubins’  TSP.  An  algorithm  for  the  Euclidean  problem  will  tend  to  schedule  very  close 
points  in  a  successive  order,  which  can  imply  long  maneuvers  for  the  aircraft.  We  will 
come  back  to  this  point  in  section  5.  In  the  following  we  suggest  a  simple  approximation 
algorithm  which  does  not  rely  on  the  Euclidean  solution. 

The  major  difficulty  is  to  determine  the  headings.  We  argue  however  that  we  have 
much  freedom  in  selecting  them,  provided  we  can  be  satished  with  a  relatively  weak 
performance  bound.  A  Erst  deterministic  algorithm  can  be  described  as  follows: 


1.  Fix  the  headings  at  all  points  to  be  0. 

2.  Compute  the  n{n  —  1)  Dubins  distances  between  all  pairs  of  points. 

3.  Construct  a  complete  graph  with  one  node  for  each  point  and  edge  weights  given 
by  the  Dubins’  distances. 

4.  We  obtain  a  directed  graph  where  the  edges  satisfy  the  triangle  inequality.  Compute 
the  solution  of  the  asymmetric  traveling  salesman  problem  on  this  graph,  using  the 
log  n-approximation  algorithm  of  Frieze  et  al.  [4]. 


The  complexity  of  the  three  hrst  steps  is  The  algorithm  for  solving  the  ATSP 

runs  in  O(n^),  so  overall  the  running  time  of  our  algorithm  is  0{n^). 

To  analyze  the  performance  guarantee,  we  use  the  bound  (1).  At  each  point,  we  can 
make  an  error  up  to  |5|  =  vr  at  each  point.  Thus  the  bound  on  the  length  becomes: 


Pij  <  (  1  +  max 


D 


(2) 


where  Dmin  =  niinj^j{(ijj},  and  C  is  defined  by  the  equation. 

Note  that  theoretically  there  is  no  reason  for  D^in  to  be  bounded  from  below,  and 
therefore  the  upper  bound  (2)  can  be  arbitrarily  bad.  In  practice  however,  in  the  context 
of  aerial  surveillance,  we  can  restrict  our  study  to  the  situation  where  the  points  have  a 
minimum  distance  between  them,  for  example  equal  to  the  coverage  radius  of  the  sensors 


of  the  aircraft,  which  allow  very  close  points  to  be  observed  at  the  same  time.  This  can 
be  justihed  by  the  following  greedy  procedure.  Suppose  we  can  observe  at  each  instant 
the  area  inside  a  disk  of  radius  Dmin  around  the  aircraft.  For  a  given  set  of  points,  start 
by  picking  a  point  arbitrarily,  and  discard  all  the  points  which  are  within  a  distance  Dmin 
of  this  hrst  point.  Next,  pick  a  second  point  arbitrarily  among  the  remaining  points,  and 
continue  similarly,  until  all  the  points  have  been  considered.  The  points  selected  have  a 
distance  greater  than  Dmin  between  them,  and  hnding  a  trajectory  visiting  these  points 
is  enough  to  cover  the  complete  initial  set.  Therefore  it  is  enough  to  run  our  algorithm 
for  a  set  of  points  with  minimum  pairwise  distance  bounded  from  below  by  Dmin,  and 
this  implies  that  for  practical  purposes  Dmin  can  be  considered  as  a  constant  independent 
of  n. 

Call  OPT  the  optimal  value  of  the  Dubins  TSP  and  a*  the  corresponding  optimal  per¬ 
mutation  specifying  the  order  of  visits.  We  have  OPT  =  P<T*{i)a*{i+i)  +Pa*{n)a*{i)  ■  = 

L{{pij},  a*),  where  the  dehnition  of  the  functional  L  should  be  clear  from  the  equation. 
Considering  the  permutation  a*  for  the  graph  problem  (where  the  edge  weights  are  the 
perturbed  distances  {pij})  and  a*  the  optimal  permutation  for  the  graph  problem,  we 
have 


Now  on  the  graph  with  weights  pij,  we  can  solve  the  traveling  salesman  problem  in 
polynomial-time  with  an  approximation  ratio  of  log  n,  so  calling  a  the  corresponding 
solution  we  have: 

L{{pij},a)  <  \ognL{{pij},a*)  <  C  log  n  L{{pij},  a*)  =  {C  logn)  OPT. 

Therefore,  we  obtain  with  the  specihc  assignment  of  headings  (or  any  assignment  in  fact) 
an  approximation  guaranteed  to  be  within  a  factor  -|-  max  y  logn  of  the 

optimum. 

We  can  make  several  remarks  about  this  algorithm: 

•  If  we  randomize  the  heading  assignments,  we  obtain  a  better  guarantee  in  expecta¬ 
tion,  more  precisely  we  can  lower  the  constant  in  front  of  the  logarithm,  by  avoiding 
the  worst-case  analysis  assuming  6  =  n.  This  is  described  in  section  4. 

•  A  deterministic  generalization  is  to  add  more  discretization  levels  to  decrease  the 
value  of  the  error  6.  At  each  point,  we  consider  a  set  of  K  possible  headings,  and 
we  want  to  select  one  of  them  optimally.  To  each  point  in  the  original  problem, 
we  associate  a  cluster  of  K  nodes  corresponding  to  the  K  different  headings.  We 
then  have  to  solve  a  problem  called  the  “generalized  asymmetric  traveling  salesman 
problem”,  i.e.  hnd  a  tour  through  n  sets  of  K  vertices,  visiting  one  node  in  each 
cluster  (i.e.  selecting  one  possible  heading  at  each  point).  The  generalized  traveling 
salesman  problem  can  be  reduced  to  the  traveling  salesman  problem  [7],  although 
in  practice  this  might  become  too  complex. 

•  As  we  have  pointed  out,  the  presence  of  Dmin  could  make  the  approximation  ratio 
arbitrarily  bad  in  theory.  It  seems  difficult  to  avoid  this  term  in  the  performance 
bound,  as  long  as  the  analysis  is  carried  independently  on  each  point-to-point  path 
as  we  did  above.  Consider  the  path  from  (0,  0,  6i  =  0)  to  (e,  0,  6j  =  0),  with  e  >  0  a 
small  number:  this  is  just  a  straight  line.  But  if  we  make  an  error  5  on  6*j,  as  e  — >■  0 


we  obtain  a  Dubins’  path  of  positive  length  since  the  initial  and  hnal  conhgurations 
are  different  and  the  aircraft  has  to  maneuver  to  change  its  heading.  Therefore  the 
ratio  of  the  Dubins’  distance  of  the  perturbed  path  to  that  of  the  original  path 
becomes  inhnite,  however  small  the  error  on  the  initial  angle  was. 

•  The  logn  factor  appears  to  be  inherently  linked  to  our  reduction  to  a  directed 
graph.  Note  however  that  discretization  and  reduction  of  a  kinodynamic  planning 
problem  to  a  directed  graph  formulation  is  standard.  The  asymmetry  is  due  exactly 
to  the  dynamic  constraints  (for  example,  making  a  U-turn  is  costly).  Working  with 
a  directed  graph  is  not  an  issue  as  long  as  the  graph  problem  is  simple  like  computing 
a  shortest  path,  but  the  situation  is  worse  for  more  complicated  problems  such  as 
the  TSP,  since  much  less  is  known  about  digraphs.  Reducing  the  logn  factor 
therefore  appears  to  be  quite  challenging  unless  a  better  approximation  algorithm 
is  devised  for  the  ATSP  or  a  method  departing  from  the  traditional  discretization 
approaches  is  used  (or  a  hardness  result  is  obtained). 


4  Randomized  Version 


We  can  modify  the  algorithm  described  in  part  3  to  obtain  an  interesting  randomized 
algorithm.  Instead  of  assigning  all  headings  to  be  0,  we  choose  the  headings  randomly 
and  independently  in  (— tt,  tt]  for  each  point. 

Consider  an  optimal  tour,  and  two  successive  points  Xj,  Xj  in  this  tour.  The  Dubins’ 
path  between  these  two  points  has  length  pij,  and  the  optimal  headings  are  9i  and  9j. 
Following  [6]  in  the  derivation  of  the  bound  (1),  we  know  that  if  we  make  an  error  of 
6  G  (— TT,  tt]  on  9i,  the  difference  in  path  length  is  bounded  by: 


Ap  <  pmax  3|5|  -|-  7r|  sin(-)|,  |5|  -|-  darccos 


sin  ^1/2 


(3) 


This  leads  to  the  inequality  (1)  by  taking  into  account  the  error  on  9j  as  well.  Now  (3) 
is  derived  for  a  change  from  9i  to  6^*  -|-  5  in  the  initial  heading.  Of  course,  we  do  not  know 
the  optimal  9^  so  the  natural  idea  is  to  choose  {9i  +  6)  uniformly  in  (— vr,  tt],  in  which  case 
the  error  6  is  distributed  uniformly  in  {—tt,  tt]  as  well.  This  implies  that  Ap  becomes  a 
random  variable  whose  expectation  is  bounded  by: 


E[Ap]  —  P  J  |3|(5|  -|-  7r|  sin  ^|,  |5|  -|-  4 arccos  —  |  sin  ^1/2^  j"  ^  <  6.79p. 


Replacing  the  corresponding  expression  in  (1),  we  obtain  as  a  hnal  upper  bound: 


E[pij]  <  (4) 

It  is  also  possible  to  rehne  the  bound  using  the  fact  that  a  CCC  path  has  length  at  least 
TTp  as  in  section  3,  but  for  our  purpose  Dmin  and  p  will  be  of  the  same  order  and  (4)  is 
then  enough. 

We  can  now  reproduce  the  analysis  of  part  3,  using  only  the  linearity  of  expectation. 
Thus,  we  see  that  we  have  a  randomized  algorithm  which,  given  a  set  of  n  points,  returns 
a  Dubins’  tour  whose  expected  length  is  within  ( 1  -f-  j  jQg  ^  Qf  optimum. 

V  Ef-fnin  } 


5  Computational  Experiments 

In  this  part  we  describe  some  computational  experiments  for  the  randomized  algorithm 
described  in  part  4.  The  algorithm  was  implemented  in  MATLAB.  For  the  computa¬ 
tion  of  the  ATSP  tours,  we  did  not  implement  Frieze’s  algorithm,  rather  we  used  LKH, 
Helsgaun’s  implementation  of  the  Lin-Kernighan  heuristic,  available  as  a  C-code  [5].  Our 
main  program  calls  this  routine  when  asked  to  solve  an  ATSP.  For  the  sake  of  comparison, 
we  also  implemented  the  alternating  algorithm,  using  LKH  for  computing  the  Euclidean 
TSP  tour  as  well.  LKH  has  excellent  performance  on  real-world  problems,  and  solves 
the  small  instances  that  we  considered  exactly  without  difficulty.  Therefore  our  imple¬ 
mentation  of  the  algorithm  returns  a  tour  of  expected  length  less  that  ( 1  -|-  |  times 

Y  ^min  J 

the  optimum  length  for  the  problems  considered,  i.e.  the  logu  factor  did  not  appear.  In 
practice,  for  a  given  instance  of  the  problem,  we  generated  10  tours  using  10  different 
sets  of  random  headings,  and  returned  the  best  tour  obtained.  The  difficulty  in  verifying 
the  performance  of  the  algorithm  is  that  it  is  hard  to  evaluate  the  true  optimum,  at  least 
when  using  a  naive  method  which  would  consist  in  computing  all  tours  for  all  possible 
headings  of  a  sufficiently  hne  grid  at  each  point. 

We  ran  simulations  for  different  sizes  of  point  sets;  the  sets  consisted  of  points  gen¬ 
erated  randomly  and  uniformly  inside  a  square  of  side  length  5.  In  all  simulations,  the 
turning  radius  of  the  vehicle  was  hxed  to  1.  In  order  to  compare  the  performance  of  the 
alternating  algorithm  and  our  algorithm,  for  n  hxed,  we  generated  10  different  sets  and 
compared  the  average  lengths  returned  by  the  two  algorithms  on  these  samples.  The 
results  are  shown  on  Fig.  2.  As  expected,  when  the  density  of  the  points  increases,  the 
performance  of  the  alternating  algorithm  decreases  with  respect  to  the  performance  of 
our  randomized  algorithm.  For  points  far  apart  however,  the  Dubins  problem  becomes 
similar  to  the  Euclidean  problem,  and  therefore  we  can  expect  the  alternating  algorithm 
to  perform  almost  optimally.  In  practice,  to  obtain  an  algorithm  which  can  handle  var¬ 
ious  point  conhgurations,  we  can  compute  tours  with  both  algorithms  and  choose  the 
best  one. 

For  sufficiently  dense  sets  of  points,  it  becomes  clear  that  the  ordering  of  the  Euclidean 
tour  is  not  optimal  in  the  case  of  the  Dubins’  TPS.  In  Fig.  1,  we  can  see  that  slightly 
modifying  the  top  loop  should  provide  a  better  tour,  if  we  modify  the  ordering  of  the 
group  of  3  points.  In  general,  as  previously  mentioned,  we  expect  that  the  Euclidean 
tour  will  schedule  close  points  successively,  which  may  result  in  long  maneuvers.  The 
alternating  algorithm  tends  to  create  numerous  loops  that  become  problematic  with  dense 
sets  of  points  (see  Fig.  3). 

We  conclude  this  section  by  commenting  on  the  tightness  of  the  performance  bounds 
obtained  in  parts  3  and  4.  For  random  sets  as  considered  above,  the  worst  case  bound, 
based  on  the  point-to-point  worst  case  performance,  is  far  from  the  actual  performance, 
typically  by  an  order  of  magnitude  or  two.  This  bound  is  very  conservative,  for  example 
the  term  1/Dmin  is  present  even  if  the  two  closest  points  are  not  scheduled  in  successive 
order.  This  calls  for  a  more  global  analysis  of  the  performance.  Note  also  that  in  our 
experiments,  we  simply  chose  the  points  uniformly  at  random,  without  imposing  any 
restriction  on  Dmin-  The  fact  that  the  square  can  be  covered  by  0{n)  boxes  of  side 
length  n“2  tells  by  a  pigeonhole  argument  that  =  0{n~^).  This  means  that 

our  approximation  ratio  is  Vt{y/n\ogn)  for  target  points  distributed  uniformly  in  the 
square.  On  the  other  hand,  Savla  et  ah  proved  in  [9]  that  the  expected  length  of 
the  optimum  Dubins’  tour  in  this  case  is  This  implies  that  our  upper  bound 


Figure  2:  Performance  of  the  alternating  algorithm  and  the  randomized  headings  algo¬ 
rithm  for  sets  of  points  of  increasing  density  in  a  square  of  side  length  5. 


(a)  Alternating  Algorithm 


(b)  Randomized  Headings 


Figure  3:  Paths  generated  by  the  alternating  algorithm  and  the  randomized  headings 
algorithm  for  a  set  of  30  points  in  a  square  of  side  length  5. 


on  the  length  retnrned  by  the  randomized  headings  algorithm  is  logn),  while 

experimental  results  show  that  we  should  expect  a  sublinear  bound  (see  Fig.  2). 

6  Conclusions 

We  presented  a  simple  algorithm  to  compute  efficiently  tours  within  sets  of  points  for 
the  Dubins’  vehicle,  which  models  as  a  hrst  approximation  the  dynamics  of  a  hxed- 
wing  aircraft.  The  Dubins’  paths  cannot  be  followed  exactly  by  a  real  vehicle  because 
instantaneous  changes  of  acceleration  are  necessary.  However,  our  algorithm  is  useful 
for  obtaining  a  hrst  approximation  of  the  trajectory  for  instance  sizes  that  could  not 
be  handled  by  more  more  precise  models.  Through  simulation,  we  demonstrated  that 
the  performance  of  the  randomized  version  of  the  algorithm  is  in  general  better  than 
for  the  alternating  algorithm,  which  is  based  on  computing  Euclidean  tours.  This  is 
supported  by  the  idea  that,  for  dense  sets  of  points,  the  Dubins’  metric  has  a  behavior 
departing  signihcantly  from  the  Euclidean  metric.  We  provided  a  performance  bound 
guarantee  for  the  algorithm,  however  typically  the  performance  obtained  is  much  better 
than  this  worst  case  bound.  Our  analysis  relied  on  a  point-to-point  result,  and  we  expect 
that  a  signihcantly  better  analysis  cannot  be  obtained  except  by  adopting  a  more  global 
approach.  For  the  randomized  algorithm,  it  would  also  be  interesting  to  obtain  a  result 
with  high  probability,  to  improve  on  our  expected  bound. 
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